Half-skyrmion picture of single hole doped Cu02 plane 
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Based on the Zhang-Rice singlet picture, it is argued that the half-skyrmion is created by the doped 
hole in the single hole doped high- Tc cuprates with Neel ordering. The spin configuration around the 
Zhang-Rice singlet, which has the form of superposition of the two different d-orbital hole spin states, 
is studied within the non-linear a model and the CP^ model. The spin configurations associated with 
each hole spin state are obtained, and we find that the superposition of these spin configuration turns 
out to be the half-skyrmion that is characterized by a half of the topological charge. The excitation 
spectrum of the half-skyrmion is obtained by making use of Lorentz invariance of the effective theory 
and is qualitatively in good agreement with angle resolved photoemission spectroscopy on the parent 
compunds. Estimated values of the parameters contained in the excitation spectrum are in good 
agreement with experimentally obtained values. The half-skyrmion theory suggests a picture for 
the difference between the hole doped compounds and the electron doped compounds. 

PACS numbers; 75.25.-|-z,74.72.-h 



In the phase diagram of high-temperature superconductors, there are a variety of phases: the Neel ordering phase, 
the spin-glass like phase, the d-wave superconducting phase, and the pseudogap phase.i This rich phase diagram is 
controlled by the hole doping concentration x and temperature. If we focus on the ground state properties, then the 
fundametal parameter is x. In order to understand the physics of high temperature superconductivity, it is necessary 
to figure out how to describe the doped holes. The goal of this approach is, of course, to find a unfied description of 
the doped holes over the whole range of x. 

In this paper, we consider the simplest case: We focus on the single hole doped system. The motivation is 
the following. Although d-wave superconductivity occurs in moderately doping region, and there is the intriguing 
pseudogap phase, it is quite hard to find a reliable description of the doped holes because of strong correlation 
effects between the holes and background spin fluctuations. In contrast, the physics of the undoped compound is 
well established compared with other phases. The system is described by the S* = 1/2 antiferromagnetic Heisenberg 
model on the square lattice and the ground state is the Neel ordered state. Considering one hole doping upon 
this well-established phase would be the first step to understand the effects of doped holes in the high-tcmpcrature 
superconductors. 

Experimentally it has been known that doped holes occupy oxygen p-orbitals in the Cu02 plane. Between an 
oxygen p-orbital hole and the nearest neighbor copper d-orbitals holes, there is strong correlations of forming a singlet 
pair. The resulting singlet state is called Zhang-Rice singlet.- Based on this picture, the t-J model was proposed. 

In terms of the t-J model, the single hole problem has been discussed extensively^ However, the focus is mainly on 
frustration effect induced by hopping of the doped hole. Not so much attention has been paid on the spin configuration 
around the Zhang-Rice singlet state. In this paper, we study the spin configuration created around the Zhang-Rice 
singlet within effective field theory approaches to the Heisenberg antifcrromagnet. 

As an effective theory for the S — 1/2 Heisenberg antifcrromagnet on the square lattice, the non-linear a model 
(NLcrM) has been studied extensively!^ 



where ps is the spin stiffness and Csw is the antiferromagnetic spin wave velocity. (Hereafter we take the unit of ?i = 1.) 
Here r is the imaginary time and (3 = (kBT)~^ with T being temperature. The unit vector n represents the staggered 
momemnt. This model is derived from the Heisenberg antifcrromagnet by applying Haldane's mapping.^ Theoretical 
formula for the antiferromagnetic correlation length (^af based on the renormalization group analysis of this model is 
in quite good agreement with experimentally obtained £,af- 

Another well-known approach is the Schwinger boson mean field theory (SBMFT).^ In this theory, the spin S* = 1/2 
is represented by 
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with the constraint Zj^Zj^ = 1. The components of a are Pauli matrices. The antiferromagnetic correlations 
are described by a mean field Aij = {zi^Zji — ZiiZji). This mean field describes also the pairing correlations of the 
Schwinger bosonsiSiS The quasiparticle excitation spectrum of the Schwinger bosons is given by 



CO, = ^X^ - (3) 

where 7^ = {sinkx ± sinfcy)/2 with the plus sign is for k^ky > and the minus si^n is for kxky < 0. The parameter 
A is originally introduced as a Lagrange multiplier to impose the constraint. In the mean field approximation, the 
Lagrange multiplier is taken to be uniform. For the ground state, A = 2JA. Bose-Einstein condensation of the 
Schwinger bosons at k = (±7r/2, ±7r/2) leads to Neel ordering>i£ 

In order to consider fluctuations about Bose-Einstein condensate of Schwinger bosons, it is convenient to use a 
gauge field description. Introducing a variable xq = CswT, the action reads S — J (fix (9^n)^ with g — Csw/Ps- In 

terms of complex fields Ca ^^d Co- with Co-Co- — the vector n is represented by n = C^C- Substituting this into 
the action, we obtain 



S= - I 

9 , 



(5^0(9^0 + (ca^C)' • (4) 



Performing a Stratonovich-Hubbard transformation at the interaction term, we obtain 
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S^- I (fxy\{d^~ia^)Ux)? ■ (5) 



a 

Note that Q is invariant under a local U(l) gauge transformation C Ce'^^^-* and C — > C^"*^'-^''- The field is a 
gauge field assocaitcd with this gauge invariance. As explicitly shown in Ref.^, the CP^ model is also derived from 
SBMFT. The relation between the fields C and C and the Schwinger boson fields (O is given by Cjo- = Zja for one 
sublattice and Cjo- = 'Zja for the other sublattice. 

In addition to these bosonic models, there is a fermionic theory. For = 1/2 spins, one can represent them by 
fermions, 



Si 



with the constraint fj^fja — 1. Based on a mean field theory, the tt-Aux phase was proposed by Affleck and 
Marston^ The excitation spectrum of the quasiparticles in the 7r-flux phase is given by 



ek = ±vy cos"^ kx + cos"^ ky, (7) 

where u is a mean field paramter of the theory. Fluctuations about this mean field state can be taken into account 
through a gauge field. For the two-dimensional S — 1/2 quantum Heisenberg antiferromagnet, the condition for the 
dynamical mass generation is satisfied^^ii^ and so the spectrum is modified as 



e/c = ± y^w^ (cos^ kx + cos^ ky) + m'^. (8) 

The mass m is associated with the ordered staggered moment .iS 

For the description of the doped hole, we shall assume that the hole forms a Zhang-Rice singlet with a copper hole 
which was proposed by Zhang and Ricc^ from the analysis of the d-p model, 

ja ia j 

~ X! '^3f-hdd]„Pic + h.c, (9) 

where ej^j^x/2 — ^j.j-y/2 = ~^j.j+x/2 — = 1- Here dj^ and pj^ are the hole creation operator at the copper 

d-orbital state and the oxygen p-orbital state, respectively. The vacuum is defined as Cu(3d)^° and 0(2p)^. Applying 
a canonical transformation and omitting unimportant terms, we obtain 



Hk^TJ -I^ + J^ d]adA PaP, , (10) 
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where Pjc, = J2{£}ej ^jiVtal'^- Constructing the Wannier wave functions,^ the Zhang-Rice singlet state is created by 
the following operator 



^ZR, = ^KT'^]i-^]i'^]T ' (11) 



where 




1 6xp (ik • (Rj — Rj')) 

1 — ^ (cos + cos ky) 



P^'o. (12) 



In this paper, we study the spin configuration around a Zhang-Rice singlet. We shall show that a spin texture 
which is called half-skyrmion characterized by a half of a topological charge is created around the Zhang-Rice singlet. 
The dispersion of the half-skyrmion is given by the same form as that of the quasiparitcle in the tt-Aux phase. The 
rest of this paper is organized as follows: In Sec^ we consider a static Zhang-Rice singlet. We argue that the spin 
configuration around such a static Zhang-Rice singlet is half-skyrmion. In Sec lIIII we construct the moving half- 
skyrmion state by making use of Lorentz invariance of the NLcrM. The lattice action is derived in Sec II VI We shall 
compare the dispersion of the half-skyrmion with the result of angle resolved photoemission spectroscopy (ARPES) on 
the parent compounds. The effective theory of the half-skyrmion is obtained by applying duality mapping in SeclVl 
Section IVTl devoted to Discussion. Finally, we summarize the results in Sec I VIII 



II. HALF-SKYRMION SOLUTION FOR STATIC ZHANG-RICE SINGLET 



In this section, we consider the spin configuration induced around a static Zhang- Rice singlet. The spin configuration 
for a moving Zhang-Rice singlet shall be discussed in Sec lIIII based on the result of this section. 

The creation operator for a Zhang-Rice singlet residing at the site j defined by Eq. suggests that the spin 
configuration around a static Zhang-Rice singlet is given by superposition of the two spin configurations: One is the 
spin configuration created by a spin up state at the d-orbital state and the other is the spin configuration created by 
a spin down state at the d-orbital state. We study these states separately, then we consider superposition of these 
states. We assume that the spin state of 0jcr does not play an important role. We shall discuss its possible effects in 
SeclvTl 

In order to study the spin configuration, we use the NLctM. Note that we cannot apply a linear response theory, 
such as a spin wave theory, to obtain the spin configuration. For the d-orbital spin state at the site j which is the 
same as that of the Neel ordered state of the parent compound before introducing the hole, the surrounding spins are 
not so much affected by the fixed spin state at the site j. However, for the d-orbital spin state which is in the opposite 
direction to that of the Neel ordered state of the parent compound before introducing the hole, the surrounding spins 
can change their directions which are beyond the range of the description of the linear response theory. 

We assume that before introducing the doped hole the spin state at the site j is up spin state and the staggered 
magnetization is in the z direction. Under this assumption, the spin up state does not change the directions of the 
neighborhood spins. If we include quantum fluctuation effects, then the staggered moments of neighboring spins 
would be enhanced because we fix the state at the site j to spin up and quantum fluctuations is quenched at this site. 
In other words, the state is in the subspace of eigenstates of the system with spin up state at the site j d-orbital. 
However, for simplicity we do not consider quantum fluctuation effects here and we restrict ourselves to a classical spin 
conflguration. Quantum fluctuation effects shall be taken into account through the renormalization of the parameters. 
Thus, the spin conflguration associated with the spin up state is n(R£) = -l-e^ for any site i. 

The spin conflguration for the spin down state at the site j d-orbital is non-trivial. We would like to obtain the 
spin conflguration satisfying the following boundary conditions: 

n(r) — > -fez (r — > oo), (13) 

and 

n(Rj) = -e,. (14) 
Since we consider the static spin conflguration, we are interested in the spin configuration that minimizes the energy, 

/dMVn)^ (15) 
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The analysis given below follows a general argument for skyrmion excitations in the NLcrM:-- We include the constraint 
|np — 1 through a Lagrange multiplier: 

E=^l d2r[(Vn)2 + A(|np-l)]. (16) 
Taking variation with respect to n, we obtain 

V2n-An = 0. (17) 
The Lagrange multiplier A is eliminated by using this equation and = 1: 

V^n-(n-V2n)n = 0. (18) 

By using the identitjii^ 

J d^r^[d^n±e^,inxd^n)f>0, (19) 
with exx = £yy = 1 and exy = ~^yx = 1, we find 

E>Anps\Ql (20) 

where 

Q = ^ / c^'i- n • {dxn X dyn), (21) 

is called the topological charge. The equality in Ea. (|20|l is satisfied if and only if 

a^n±e^,(nx9^n) = 0. (22) 

The solution of Ea. (|18|l is divided into sectors with different Q values. Since the energy in each sector has the lower 
bound determined by Ea. (|2()|l . it is enough to solve Ea. (|22|l for our purpose. To solve this equation, it is convenient 
to use a variable w — {ux + iny)/{\ — riz). In terms of Eg. 1)22(1 is 

dxW — —idyW, dxW = idyW. (23) 

Introducing z = x + iy and z = x — iy, we find 

d^w = 0, = 0. (24) 



Therefore, the solution satisfies the Cauchy-Rieman eqation. Analytic function of z or z is the solution of Ea. (|22|l . In 
terms of w and w, the vector n is described by 

+ w . w — w |wp — 1 



|w|2 + l' Vp + l'kP + 1/' ^^^^ 

The boundary conditions of (|13|l and H14|) are 

\w\ oo (r ^ oo), (26) 

and 

w = (r = R,), (27) 
respectively. The solution that satisfies these boundary condition is 

w^j, (28) 



or 



w = j, (29) 
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up to a phase factor that is associated with a global rotation of all spins in the plane. Here A is a parameter which 
is associated with the size of the spin configuration and Rj is taken to be the origin to simplify the expressions. The 
vector n representation of Eqs.ljSHl ^^'^ is the following 



2Xx 2Xy - 



and 



r2 + A2 ' r2 + A2 ' r2 + A2 
2Xx 2Xy r2 - A' 



The topological charge Q is given by 



r2 + A2 ' r2 + A2 ' r2 + A2 



(1+|U.|2)2 



(30) 



(31) 



Q - - / d^r^A^l^^^^A^;^^^,^^^^^. (32) 



For the solutions Eas. (|28() and H29|) . the topological charges are Q — 1 and Q = —1, respectively. These solutions are 
called skyrmion. The energy of these skyrmion solution is E' = '^'^Ps^ which is calculated from the following expression: 



j2 {dw){dw) + {dw){dw 

(1 + |W|2)2 

Note that there is no solution in the Q = sector. In fact, solutions in this sector satisfy dui = and doj — 0, and we 
obtain lu — const. Obviously such a solution does not satisfy the boundary conditions. 

Now we consider superposition of the uniform state and the skyrmion state. Unfortunately, the classical super- 
position of the two spin configuration is not the solution of the field equation. However, these solutions suggest 
that the resultant spin configuration is chracterized by a topological charge Q with < \Q\ < 1. The value of Q 
is determined by making use of the fact that Neel ordering state is described by Bose-Einstein condensation of the 
Schwinger bosons. W In order to examine the value of Q, we use the following representation of Q by the CP^ gauge 
field a a. 



I (Pr 

Q = I i^x^y - dya^) . (34) 



From this expression, we see that the spin configuration with Q corresponds to the flux 2ttQ in the condensate 
of the Schwinger bosons. Since the spin 1/2 bosons C,a{x) are confined in the Neel state, all bosons are paired in 
the low-energy physics. Because pairs of the bosons carry the gauge charge two, the flux quantum is tt similar to 
the conventional BCS superconductors. Therefore, the flux value is not arbitrary and Q must be in the form of 
Q — nn, with n being an integer. Meanwhile, from the constraint < IQj < 1, the flux associated with the spin 
configuration satisfies < 27r|(5| < 2tt. Thus, we conclude 27r|(5| = tt, or \Q\ = 1/2.^™ Since the topological charge is 
one-half, we call this spin configuration as half-skyrmion. The energy associated with the half-skyrmion spin texture 
is 47r/9s|Q| = 2ttPs = Eq because we can apply Ea. H22|) outside the core region. Due to the limitation of the effective 
theories, we cannot determine the core energy. It would be determined from a calculation based on a microscopic 
model. The half-skyrmion solution is also discussed in the ferromagnetic Heisenberg model^^ in the context of quantum 
Hall systems. Since the value of \Q\ = 1/2 is obtained from the calculation of the topological charge using Eqs. H3U|) 
and (|31() with excluding the core region r < A, we may use Eqs. (|30|) and (|31|1 for the expressions of the half-skyrmion 
and the anti- half-skyrmion spin textures, respectively. The half-skyrmion and anti- half-skyrmion spin textures are 
shown schematically in Figs. ^ and |21 

III. MOVING HALF-SKYRMION SOLUTION 

In the last section, we have argued that the spin configuration around a static Zhang-Rice singlet is the half- 
skyrmion. In this section, we construct the moving half-skyrmion solution from the static half-skyrmion solution by 
making use of the Lorentz invariance of the NLcrM and the CP^ model. 

The action |^ is written in the Euclidean space time. The form in Minkowskii space time reads 



1 / dn 



L, \ dt 



2 



(Vn)2 



(35) 



FIG. 1: The half- sky rmion spin texture. Arrows indicate the staggered moment and the filled circle at the center represents 
the Zhang-Rice singlet formed site. 
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FIG. 2: The anti-half-skyrmion spin texture. Arrows indicate the staggered moment and the filled circle at the center represents 
the Zhang-Rice singlet formed site. 



Apparently the action is invariant under Lorentz transformations with Cgw being the speed of "light." Let us consider 
a following Lorentz transformation, 

^ - (36) 
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It is easy to check that the action is invariant under this Lorentz transformation. 

Now we apply the Lorentz transformation H36fl and (|37|l to the static half-skyrmion solution (|30|) : 



(37) 



(38) 



The excitation spectrum of the half-skyrmion is obtained by calculating the energy-momentum tensors. We write 
Ea. H38|l as n' = ns(x', y), where x' = ^{x — vt) with 7 = 1/^1 — (u/csw)^- The energy-momentum tensors are given 

by 

Tj; = psd^nd^n - ^d^ndpuSi^, (39) 
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where A^^Bn = AqBq — A ■ The energy is 



E = / d^rT! 



2 



^ / 



4^ + 7M (a,,n,)V(a,n,)2 



7Ps j dxdy{dx'nsf 



(40) 



where we have used that dtn' — —^vd^'n-g, d^n' = "/d^'n-s, dx = jdx' , and (d^'ng) — (dyiig) . The x-component of 
the momentum is calculated as follows 



P. 



± I d^rT? 

Ps 



dxdy (dtUs) ■ (S^n^) 



(41) 



Pn = 



d\T^ 



dxdy{dtns) ■ {dyiis) 



From Eqs.lgni), l|lT)l . and we find 

E^=cl^P^+El 

By considering general Lorentz transformations, we find that the following relation holds: 

E'^c',^{P^ + P^)+El 
Therefore, the excitation spectrum of the half-skyrmion is given by 



Ei^ 



± 



El 



(42) 
(43) 
(44) 

(45) 



Having obtained the relativistic dispersion (|45|l . we consider the action of the half-skyrmion. Since the half-skyrmion 
is a topological spin texture, one might expect that Berry phases affect the statistice of the half-skyrmion as in the 
fractional quantum Hall systems. In the fractional quantum Hall systems, the effective theories are characterized by 
topological field theory. Berry phase effect determines the statistics of quasiparticles. In constrast, the system of 
single hole doepd antiferromagnet is not characterized by topological field theory. Indeed, a gauge field that describes 
the Berry phase effects is massive due to Bose-Einstein condensation. Therefore, the leading term of the gauge field 
is the mass term. In such a situation, we do not expect that Berry phases play an important role in determining the 
statistics of quasiparticles. In the absence of the Berry phase effects, the statistice of the half-skyrmion is fermion 
simply because the doped hole that obeys the fermionic statistics sits at the core. 

The statistics of the Zhang-Rice singlet is infered from its field operator. Since the d-orbital hole states constitute 
localized spin 1/2 moments and one can choose either a fermionic description or a bosonic description for the spins, 
the statistics of the Zhang-Rice singlet is either fermion or boson. However, now we are interested in the Neel ordering 
phase. Bosonic descriptions, such as Schwinger boson mean field theory or NLctM, are suitable for the description of 
the Neel ordered state. In this case, the statistics of the Zhang-Rice singlet is fermion. 

Fermion field obeying the relativistic excitation spectrum (|45() is described by a Dirac fermion action. The action 
of the half-skymion may be written as 



C = 



s=± 



(46) 
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with ipg = tpl^o and mCg^ = Eq. The index s is for the sign of the topological charge. That is, s = + is for 
the half-skyrmion and s = — is for the anti-half-skyrmion. The Dirac fermion ip has four components: There are 
positive and negative energy states. These states obey the dispersion (|45(l with the origin either at ki = (7r/2,7r/2) 
or k2 — (— 7r/2, 7r/2). This is suggested from the fact that the Schwinger bosons are gapless at four zone centers 
(±7r/2, ±7r/2). Because Q ~ (tt, tt) connects two points in the diagonal directions, there are two independent points 
ki nad k2. The 7 matrices are 4x4 matrices and satisfy 7^71/ + 71^7^ = 2(5^,^. 



IV. HALF-SKYRMION ON THE LATTICE 



The expression on the lattice is derived by discretizing the continuum Lagrangian Ea.H46|l. After Fourier transfor- 
mations, we obtain 



"^cL + dr sin +i sin ky 1 



sin k^ + i sin k,. 



dr 



a=l,2 a 

By shifting the origin in the momentum space to either ki or k2, we obtain 

mCg^ + dr cos fcj, + i cos k, 

COS k,j: + i COS k 



(47) 



^-E 



V'i.(fc) 



sw 

■ COS kx -\- i COS ky 



1pla{k). 



mc^„ + dr 
COS k^ + i COS ky 



dr 



dr 



Diagonaliztion of the matrix leads to the following expression, 



±\/ Csw(cos2 kx + cos2 ky) + {mcl^Y 



(48) 



(49) 



Note that there is no fermion doubling problem, which occurs when one formulates a Dirac fermion on a latticed 

The parameters Cg^ and mc^^ = 2'Kps are determined from the values for the Heisenberg antiferromagnet. We 
use Zc = 1.17 and Zp = 0.72 which are estimated from quantum Monte Carlo simulations^^ and a series expansion 
technique. Substituting these values into Eq. H49|l . we find that the band width is ^ 1.5J and m/ J ~ 1.13. Meanwhile, 
experimentally estimated band width by Wells et al. is ^ This discrepancy would be associated with the 

deviation of the real system from the NLctM. For the mass value, Ronning et al. evaluated it by using the form 
— mc^^. The result is m/ J ~ 1.3i^ Some of this discrepancy might be associated with the mass renormalization 
due to antiferromagnetic spin fluctuations. This shall be discussed in SeclVl 



V. EFFECTIVE THEORY FOR HALF-SKYRMION 



In order to include antiferromagnetic spin fluctuation effects on the half-skyrmion, we shall derive the effective 
theory of the half-skyrmion. We make use of a duality mappingS^ for that purpose. 

Before the application of the duality mapping, we point out that a half-skyrmion can be seen as a vortex in the 
CP^ model. As stated in Sec^l the topological charge of the half-skyrmion is represented by the gauge flux with 
respect to the CP^ gauge field a^. In the CP^ model, the vector n reads 

n = (CtCi + CiCT,-* (CtCi - CiCr) ,CtCt - C|Ci) • (50) 

Since n is a unit vector, C and C, satisfy C|Ct + CxCi = 1- In terms of Co- and Co-: the half-skyrmion solution Eas. H3Ufl 
and H31|l has the following form 

•^T = ; / exp(±z6>)exp(zx), (51) 
Vf + 

Ci - ^=i=cxp(zx), (52) 

where the minus sign is for the half-skyrmion and the plus sign is for the anti-half-skyrmion and % is a constant. In 
general, the half-skyrmion solution is represented by 



Ci ) \v u* )\ (r/Vr2 + A2) exp {±i6 



(53) 
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where u and v are constant complex numbers and satisfy |u| 



= 1. The matrix 



is a global SU(2) 



transformation. The boundary condition at infinity is transformed to n ^ (— mw* — vu* , ~i[uv* — vu*), — |up + |tip) . 

In order to see the relation between the half-skyrmion and a vortex, we take u ~ v = l/\/2- In this case, the 
half-skyrmion solution have the following form 



Ct 



A — r exp {±iO) 
V2(r2 + A2)"' 

A + r exp {±.i9) 
v/2(r2 + A2)"' 



This has a vortex form at r ^ A: 



exp(±j6 



l/%/2 



(54) 
(55) 

(56) 



Thus, the half-skyrmion is seen as a vortex introduced in the system. Such a vortex is taken into account in the CP^ 
model as follows 

Ca = pI'^ exp(z0) = pI'^ exp(z0o + i(l)v)- (57) 
Here describes coherent motion of the bosons and 0„ describes the vortex. In case of Ea. ()56|l . (j)^ is 



±tan"^ 



(58) 



with {xv,yy) representing the coordinate of the vortex. 

We rewrite the CP^ model by substituting Ea. (|57|l into Eq.jSJ: 



(59) 



The amplitude fluctuations would be important only in the vicinity of the core. We focus on the outside of the core 
and assume a constant value for p„: po = TliaiPi^)- After introducing a Stratonovich-Hubbard field J^, we obtain 



4po 



(60) 



The field is associated with the spin current. Integrating out 0o leads to d^J^ = 0. From this equation, we can 
represent in terms of a gauge field 



Zn 

Substituting this equation into the action, and integrating out a^, we obtain 

Sd = Sa + Sint, 

where 



1 

4^ 



and 



with 



/ d^^A^Jl 



(61) 



(62) 



(63) 



(64) 



(65) 
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The value of the gauge charge ca depends on the action of the gauge field , which would be massive because of 
Bose-Einstein condensation of the Schwinger bosons. Here we treat ca as a parameter of the theory. 

Since describes the spin current and A^j, is related to through Eq. l|HT|l , the gauge field is associated with 
spin excitations, such as antiferromagnetic spin waves. Indeed, is massless in the Neel ordered phase and the 
velocity of the mode is Cgw 

By taking into account the fact that the half-skyrmion is given by Ea. (|46|l . we obtain the following action for the 
half-skyrmion 



S ^ J tfx ^i)^ [7^ {d^ - iq^A^) + m] ipa + {df^A^ - d^A^) 



(66) 



where qa is an index for the sign of the topological charge. Hereafter, we set Csw = 1- 

In order to study the effect of spin fluctuations on the half-skyrmion, we formulate the theory on the square lattice. 
The lattice form of H66|l is given by 



a j 



3 



We expand exp(±igo-0^(i)) with respect to <f)^{i). The first order term C^^^ is 

-^^'^ = E E E l^'^k+.^alMD^Kae-^"'"'^ cos (k^ + 1^ 

a ^ k,q 



The second order term is 



^^'^ - -^E EE W,Va7.'^.(9)^.(9')V'..ae-^(^^+«^'^/^ Bin ffc, + 



k.q.q' M 



(67) 



(68) 



(69) 



We evaluate the self-energy S^^^ that comes from the second order term of C^^^ and S^^^-* that comes from the first 



order term of £^^^ . The self-energy ' is given by 



^^^^^ " 2T]E Dg-f^aink^ 



(70) 



with Dq — 1/ [{iujn)^ — LUg] . The mass renormalization due to this self-eneryg is 



m 



l-0.123eV 



The self-energy S^^-* is given by 



yil) _ 2 



-i^DqGk+q^^i cos^ {k^_, + 



(71) 



(72) 



Substituting the explicit form of Gfc, 



we obtain 



'^k = (-«7m sin k^ + m) . 

sm Ku + TO^ 



(/) 2e\ [ijp sin(fcp + gp) + mc^] cos^ (fc^ + ^) ~ 2^7At sin(fc^ + ^At) ' 

2.,, - 



d?q 



sin^ qu [sin^(fci. -|- q^) + rn?] 

EpCOS^ (fc^-(-f) 



(2^^)' (E, sin2 q,) ^ sin2(fcp + q^) + 



(73) 



(74) 
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Since the dominant contribution comes from the region of and 5^ (t'',''!'), we make an approximation for 

evaluating the integral: 



p 

Em cos^ fc^ 



TT sin^ fci. 

At the zone centers, the energy is shifted by —2e\/{-Km). 



(75) 



VI. DISCUSSION 



On the description of the half-skyrmion spin texture formation, it is essential that the Zhang-Rice singlet state 
has the form of superposition of the d-orbital spin up state and the d-orbital spin down state. However, each of 
d-orbital spin states accompanies the oxygen p-orbital hole state. As argued by Aharony et al.,^^ there is a possibility 
of replacing the antiferromagnetic superexchange interaction with a ferromagnetic interaction on the bond where the 
doped hole occupies the oxygen p-orbital state. In the picture of the Zhang-Rice singlet where a copper site hole spin 
forms a singlet state with the Wannier state of the symmetric combination of the four oxygen p-orbital hole states 
surrounding the copper site, the interactions between the nearest neighbor spins could be replaced with ferromagnetic 
interactions. If this is the case, the skyrmion spin texture would be formed for the opposite spin state compared with 
the antiferromagnetic interaction case. However, obviously this does not affect the conclusion on the half-skyrmion 
formation. Another possibility on the effect of the oxygen hole state would be related to the core structure of the 
half-skyrmion. The discussion on the core structure is beyond the description of the effective theories. For the 
investigation of this effect, one would require a microscopic model, such as the d-p model. 

As argued in Sec0 the half-skyrmion spin texture can be seen as a vortex introduced in Bose-Einstein condensate of 
the Schwinger bosons that describes the Neel ordering. Generally, vortices disturb the phase coherence of condensate. 
Rapid suppression of Neel ordering by hole doping could be understood by this pictured We expect that this effective 
Neel order suppression is absent in the electron doping case. Our picture for the half-skyrmion formation is based on 
the fact that the Zhang- Rice singlet has a form of superposition of the spin up and spin down states. This suggests 
that a picture for the difference between the hole doping and the electron doping. In the electron doping case, the 
doped electrons occupy copper sites. The occupied copper site has spin zero. However, the electronic state at this site 
does not have a form of superposition of different spin states but is simply given by a doubly occupied site. Therefore, 
we do not expect the half-skyrmion formation for the electron doped systems. Similar arguments can be applied to 
Zn doping. If we replace the copper with Zn, the electronic state at that site does not have superposition of different 
spin states. Since vortices are not induced by either the electron doping or the Zn doping, those dopings are not so 
effective to destroy the Neel ordering as observed experimentally. By constrast, if we replace a copper with Li, then 
a hole is introduced. That hole is believed to occupy an oxygen p-orbital state. If we assume that the Zhang-Rice 
singlet formation between the hole spin and a copper site spin, the half-skyrmion would be created. This is consistent 
with the experiments that report that the critical Li doping concentration for the destruction of the Neel ordering is 
almost the same as that of the hole doping concentration.^^ A skyrmion-like spin texture^^ formation for the Li-doping 
case was discussed by Haas et alM 

In Sec. IIVI we have argued that the dispersion H49() is qualitatively in good agreement with the ARPES result 
on the parent compounds. A similarity between the dispersion 149|) without the mass term and the ARPES result 
was first pointed out by Laughlin.'^i However, in the absence of the mass term cusp structures appear around the 
(±7r/2, ±7r/2) points which does not agree with the experiment. This discrepancy is corrected by including the mass 
term. Recent estimation of the mass termS reports a value which is close to the theoretical prediction as discussed in 
Sec. IIVI Another point of view has been suggested based on a self-consistent Born approximation- of the t- J model. 
The analysis of the t-J model predicts a relatively flat dispersion along the (tt, 0) point to the (0, tt) poiint whose band 
width is much smaller than the ARPES result. This discrepancy is improved by taking into account next and third 
nearest neighbor hopping terms. Contrastingly, the dispersion (|49() predicts the same dispersion along the (0, 0) point 
to the (tt, tt) point and the (tt, 0) point to the (0, tt) point. This seems consistent with the experiment by Ronning et 
alM. 

The quite broad peaks observed by ARPES would be associated with the coupling to some bosonic modes. The 
question is what boson mode plays the major role for this broadening. A scenario based on an electron-phonon 
coupling has been proposed by Mischenko and Nagaosai^ It is argued that the electron-phonon coupling is in the 
strong coupling regime for the quasiparticles in the t-J model. The quasiparticle dispersion in the t-J model is 
associated with the antiferromagnetic spin wave effects. Whereas in the half-skyrmion picture, the dispersion is given 
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by the soliton character of the half-skyrmion in the spin system. The antiferromagnetic spin fluctuations can play 
some role for the broadening of the spectrum. This issue will be examined in the future publication. This effect 
comes from the coupling between the half-skyrmion and the gauge field that is associated with the antiferromagnetic 
spin fluctuations as discussed in Sec. |V| Of course, this just suggests another possibility for the broadening of the 
spectrum. 

Ng argued^ that the vortex excitations in the Schwinger boson mean fleld theory correspond to the quasiparticles 
in the 7r-flux phase. Similarities between the skyrmion-like spin texture^ and the quasiparticle of the 7r-flux phase 
was pointed out by GoodingS^based on a numerical simulation. These seem to be consistent with the half-skyrmion 
picture. In the context of the spin-charge separatiouj^^ Baskaran argued^! that the half-skyrmion can be seen as a 
deconfined spinon from the analysis of the n-skyrmion solution of the NLcrM. 

The half-skyrmion picture may be extended to the slightly doped regime. If we increase the number of holes, then 
the interaction between the half-skyrmions becomes important. A half-skyrmion lattice state might be possible in an 
appropriate doping concentration. We expect that a skyrmion-like spin texture is formed even in the magnetically 
disordered phase. Since the antiferromagnetic correlation length ^af is finite, the topological charge that characterizes 
the spin texture is given by 



Probably the term "topological" is not appropriate in the disordered regime because the index Q would be no longer 
quantized. If doped holes accompany a spin-texture with non-zero Q, then we can apply a mechanism of d-wave 
superconductivity based on a skyrmion like spin texturei^? 



To summarize, it has been argued that the half-skyrmion spin texture is created by doping a hole into the Cu02 
plane. The picture is based on the Zhang-Rice singlet. The Zhang-Rice singlet wave function has the form of 
superposition of different d-orbital spin states. The half-skyrmion formation is the result of superposition of the 
skyrmion spin texture and the trivial uniform state associated with each of d-orbital spin states. 

The excitation spectrum of the half-skyrmion is qualitatively in good agreement with ARPES experiments on the 
parent compounds. Although the theory is based on the effective theories, such as the NLctM and the CPi model, 
estimation of the parameters in the dispersion gives us the values close to experimentally obtained values. 

Since the half-skyrmion picture is based on the superposition nature of the Zhang- Rice singlet wave function, the 
picture is not applicable to the electron doped systems where the description of the doped electrons is not expected 
to have that character. 

Interesting extention of the half-skyrmion picture would be to examine whether the picture can be applicable to 
slightly doped regime where the interaction between the half-skyrmions is not negligible. It would be also interesting 
to examine whether skyrmion-like spin texture is formed in the disordered spin state. In the absence of the Neel 
ordering, the topological charge could be arbitrary. However, we expect that strong antiferromagnetic correlations 
can stabilize a skyrmion-like spin texture created by a doped hole with non-vanishing topological charge. 
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